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Abstract 

Quantum vortex tunneling is studied for the case where the Hall and the 
dissipative dynamics are simultaneously present. For a given temperature, 
the magnetization relaxation rate is calculated as a function of the external 
current and the quasiparticle scattering time. The relaxation rate is solved 
analytically at zero temperature and obtained numerically at finite tempera- 
tures by the variational method. In the moderately clean samples, we have 
found that a minimum in the relaxation rate exists at zero temperature, which 
tends to disappear with increase in the temperature. 
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I. INTRODUCTION 



The depinning properties of vortices in high temperature superconductors(HTSC) have 
generated a good deal of interest over the past decade. [|l],^| Yeshurun and Malozemoff || 
reported on the existence of the giant flux creep which arises from the thermally activated 
motion of vortices from one metastable state to a neighboring one. The probability for 
such a process is proportional to exp(— [T/fceT), where U is the height of the energy barrier 
which depends on the pinning strength and the external current. [Q At an extremely low 
temperature the exponent diverges and the vortex cannot move out of the pinning sites 
any more. Hence, the dynamical magnetization relaxation rate Q defined as ksT/U is 
expected to vanish at T = 0. However, many experiments [|5HlO[l have demonstrated that 
the relaxation rate does not disappear at sufficiently low temperatures, which leads to the 
existence of quantum tunneling of vortices trapped in the pinning potential. 

In general, quantum vortex creep is well described by the dynamics of two major forces: 
the Hall force and dissipative force. Within the collective pinning theory, Blatter et al. 



11] considered the quantum vortex tunneling for the case where the dissipative term is 



dominant in the motion of vortices. On the other hand, Feigel'man et al. [|12j proposed that 
the Hall tunneling is dominant in clean superconductors by estimating the low-lying level 
spacing in the vortex core and the transport relaxation time of the charge carriers. Many 
experimental results have been interpreted within the two frameworks. Recently, however, 
van Dalen et al. |§ observed experimentally that the vortex tunneling in HTSC may occur 
in the intermediate regime between the purely dissipative tunneling and the superclean Hall 



tunneling. Feigel'man et al. []12] and Morais Smith et al. |y| studied the problem in the two 
regimes, but they only obtained the qualitative results based on the scaling analysis of the 
action. The main difficulty of the problem is in the fact that there is the time nonlocality 
caused by the dissipative dynamics. Recently, the present authors |14| have treated the 
problem quantitatively by using the variational method and presented the numerical results 
for the magnetic relaxation rate at zero temperature in the intermediate regime. Later, 
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Melikidze |T5) studied a similar problem by considering the quadratic Hamiltonian of the 



vortex coupled to the environment. Through the analytic diagonalization, he obtained 
the dynamical magnetization relaxation rate at zero temperature as a function of the Hall 
and dissipative coefficients and found the minimum feature in the intermediate regime. 
Previous works have treated the problem only at zero temperature, but we extend it to 
finite temperature in this work. Based on the instanton approach, we have obtained the 
numerical results for the relaxation rate at finite temperatures and its analytic expression 
around the crossover temperature between thermal activation and quantum tunneling. Using 
the functional dependence of the relaxation rate on the Hall and dissipative coefficients at 
the crossover temperature, we have also obtained the analytic expression of the relaxation 
rate at zero temperature. 

This paper is organized as follows. In Sec. II, we introduce the general formulation for the 
vortex tunneling rate in the presence of the Hall and the dissipative dynamics based on the 
instanton method, and discuss the Ohmic dissipation formulated by Caldeira and Leggett 
]T6| . In Sec. Ill, we calculate the magnetic relaxation rate by taking into account the pinning 
potential barrier generated by impurities. Writing the action and the corresponding classical 
equations in the Fourier space, we analytically calculate the relaxation rate as a function 
of the external current and the Hall and the dissipative coefficients at zero temperature. 
We also discuss the minimum of the relaxation rate in the intermediate regime. In Sec. 
IV, we numerically calculate the finite-temperature relaxation rate based on the variational 
method. It is found that the minimum in the relaxation rate tends to disappear with the 
increase in the temperature. We conclude in Sec. V. 

II. BASIC FORMULATION 

We consider the pancake vortex in the xy plane with length L c along the z axis. L c is the 
collective pinning length which can be expressed in terms of the mass anisotropy parameter 
e \ = m /M < 1, the coherence length £, the depairing curremt density j , and the critical 
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current density j c : L c ~ £ a £,(jo/ jc) 1 ^ 2 , within the weak collective pinning theory. |IJ] L c is 
obtained by minimizing the energy density which includes the elastic energy of the vortex 
string, the energy gain from the random pinning potential and the contribution from the 
Lorentz force. Thus, each segment of the length L c of the vortex is pinned by the collective 
action of all the defects within the collective pinning volume V c ~ ^ 2 L C . 

To study the quantum tunneling of the pancake vortex at a finite temperature, we 
consider the path integral representation of the partition function given by 



Z((3h)= <f D[u(r)]ex P (-S E /h), 



(1) 



where (3 = 1/kgT and Se is the Euclidean action. The path sum includes all the periodic 
paths u(t) = u(r + (3fi), where u is the displacement vector of the vortex in the xy plane. 
The Euclidean action Se includes the Euclidean version of the Lagrangian Le : 



S e [u(t)} = / drL E [u(r)]. 
Jo 



(2) 



The tunneling rate F in the semiclassical limit, with an exponential accuracy, is given by 



T ocexp[-S™ n (T)/ft]. 



(3) 



We study S^ m (T) which gives the trajectory with the period (3h that minimizes the Eu- 
clidean action. jl8| Considering the situation where the inertia term is not relevant and the 
vortex dynamics is dominated by the Hall and the dissipative forces, we write the Euclidean 
action as 



/ dr{L 
Jo 

+E 



—ia—^-Uy + V(u x , Uy) 



I i 

-m k {xl + yl) + -m k u 2 k 
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+ Vk 



m k u z k , 



(4) 



where a is the Hall coefficient and V(u x , u y ) is the pinning potential per unit length which 
includes the contribution from the Lorentz force. The last term of Eq. (fj) represents the 
dissipative environment of the vortex consisting of a set of harmonic oscillators as formulated 



by Caldeira and Leggett. [16] The effect of the dissipative environment is characterized by 
the spectral function 



z k mkiOk 

With the oscillators integrated out, the Euclidean action takes the form 



(5) 



dr{L c 



+ - 



dr'K (r - r') \{u x {r) - u,(r')) 2 + (u y (r) - S (r')) 2 l}, 



(6) 



where the nonlocal influence function is expressed as 



k o(t) = — / duJ(u)exp(-u)\T\). 
In 



(7) 



III. QUANTUM TUNNELING OF A VORTEX 



In order to study the motion of a vortex, we need to first analyze the structure of the 
model potential V(u x , u y ). Since the external current j along the y direction brings the 
system into a metastable state by tilting the potential, the vortex has a chance to move out 
of the pinning potential. Let us define u Xi as the critical position of the vortex at which the 
barrier vanishes at the critical current j c . In the limit j — > j c , u Xi and j c satisfy 



~dV~ 




'd 2 V~ 


du x 


u x =u x . 


du 2 x 



0. 



(8) 



With V{u Xl u y ) = Vp(u x ,u y ) — 4>oju x /c, u Xi and j c are given by the relations 
(dV p /du x ) Ux=Uxi = 4>oj c /c and (d 2 V p /dul) Ux=Uxj = 0, where O — hc/2e is the flux quantum. 
For the pinning potential, we choose an appropriate model potential describing a typical 
tunneling situation: V p (u x , u y ) should exhibit a local minimum and should be connected via 
a saddle to the free space along one direction (we choose this direction as the x axis.). A 
model potential satisfying the requirement is 



V(U x ,Uy 



2 



u , 



R 



RJ 



u. 



R 



(9) 



where c x = R 2 M^(cV^\(d 3 V/dul) Ux=Uxi \V 2 , c 2 = R 3 \{&V/dul) Ux=u \/{2V ), and e 



1 ~~ j/jc <S 1. In Eq. (H), Vo and R are the height and the range of the pinning potential, 
respectively, and for a typical weak pinning potential Vo — (4>o/4:Tr\ X y) 2 and R(> £). And 
ci_2 are the dimensionless coefficients of the order of 1 and \ xy is the bulk-planar penetration 
depth. 

To consider the tunneling of a vortex in the two regimes, we investigate the behavior of 
the Euclidean action (Q) . In order to estimate the order of magnitude of each term in the 
action and to simplify the calculation for e < 1, we introduce the dimensionless variables 

( 2 ° 2 \ - ( 2 ° 2 \ - ( V^ V o \ /t n \ 

«. = (^gK. °. = (^5^K T = (7f^) T - < 10 > 

where «o — 7r7in a and n s is the number density of the electrons in the condensate. 

Assuming the Ohmic dissipation where the frictional force acting on the vortex is 
linear to the vortex velocity, fll Ul9| , |20| the spectral density becomes J(oS) = rjuj, where 
r\ = (7r/2)J2i{Cf/miUjj)5(u — ivi) =constant. |L6| With this choice, we have the influence 
function 

A'oW = (ID 



which leads to the Euclidean action 



S E = (^^)(L c a R 2 y/ 2 I RD , (12) 



5/2 



where 



,A 

-^hd = / dr{—ia 
Jo 



du x _ 1 _ 2 1-2 1 



4 J-oo \f — fi\ 2 

where A = /3frVo v /eci/(v / 2ao-R 2 )- The dimensionless Hall (= a/(\/2a )) and dissipation 



coefficients (= rj / \L c \/2 : K 2 c^eaQ)) are given by pO 



V^l + ^or,) 2 ' 1 ' 

Vl ~ 2n v ^l + (o; r r ) 2 ' 1 J 



Here, uo is the level spacing of the quasiparticle bound states inside the vortex core and 
r r (= m/ne 2 p n ) is a quasiparticle scattering time, where n is the number density of the charge 
carriers and m and p n are their effective mass and resistivity, respectively. As can be seen 
in Eqs. (fL4j) and (|15D, the Hall coefficient a is reduced from its pure value 7rhn s due to the 
dissipative effect. Although Ao et al. suggested that the Hall coefficient is originated from 



the topological property and thus not renomalized, it seems that at least some aspects 
of the experimental behavior || can be understood on the basis of the renomalization of the 
Hall coefficient. Therefore, it is meaningful to take a and 77 to be two parameters determined 
by the magnitude of u T r . 



A. Action in the Fourier Space 

When the Hall and the dissipative dynamics are simultaneously present, the classical 

trajectories of u x and u y satisfy 

du v ui r°° , ,du x . 1 . , 

* a 'iF + ^-f-^L^<^-?= ' (16) 

du x _ f°° ,du v , 1 , . 

-%a x — + u y - 7]ic x e / dri{— -)- r = 0. (17) 

CLT J— 00 &T\ T\ — T 

The substitution f — > — f in Eqs. (|16]) and ( |T7| ) shows the invariance of the equations by 
taking u x (—f) = u x (f) and u y (—f) = —u y (f). We will keep c\ in the ensuing equations, 
although we will take c\ = 1 for the numerical calculations. Denoting u(r) = (u x (f),u y (f)), 
we have u(f + A) = u(f) at a finite temperature. A simple analysis shows that u x (f) is real 

and u y (f) pure imaginary, so they can be expanded into the Fourier series: 

00 

Ux(t)= u n exp(iu) n f), (18) 



n=— 00 
00 



(r) = -% v n exp(iuj n T), (19) 



where u n = 2nn/A (n=0, 1, 2...). Substituting them into Eqs. (|T6|) and QT?D, we have 

( ,_ , a l^l \ 1 ^ 

1 + 7T?7i \UJ n \ + — — — \u n = - > u n+m u m , (20) 

V 1 + 7re77iCi|^ n | / 2 ' 



m=— 00 



where U- n = u n and t>_ n = —v n . Although Eq. ( |20| ) is a one-dimensional problem with 
respect to u n , its solution becomes complicated by the presence of the nonlocal term arising 
from the cubic potential. For general ujqt t and e, we have numerically solved Eqs. ( |20| ) 
and ( pip via the variational method. The trial function for the variational method has 
been taken by combination of the analytic solutions in the two extreme limits as follows. 
At zero temperature, u x (uys are of the form Co/ sinh(o;i7ra)) in the Hall limit {uo$T r — > oo) 
and exp(— 71771 |u)|) in the dissipative limit {ujqt t — > 0), 0] so a natural choice for the trial 
function at finite temperatures is 

M « = ■ rr~^ x + P3exp(-p 4 \u! n \), (22) 
smh (p 2 uj n ) 

where p^'s (i=l,2,3,4) are free parameters to be determined by the variational method. It 
turns out that the numerical variational method with the trial function works very success- 
fully. 

Using the Fourier series in Eqs. ([T5|) and ([19|), we write /hd as 



'HD 



«l^ n M n f n + - ~vl - ^U n ( E U n+m U m \ + ^?7i|^ n |(^ - Ci€V 2 n ) 



\m=—oc 



which further reduces to 



(23) 



1 00 / n/V; 2 \ 

J hd = -A £ 1 + vr^l^l + — <- (24) 

6 „,±^oo\ 1 + 7re?7iCiK|/ 



B. Quantum Relaxation near the Crossover Temperature and at Zero Temperature 

At T c , the crossover temperature between thermal activation and quantum tunneling, 
the classical trajectories become independent of f, i.e., u x (f) = 2 and u y (f) = 0. When 
A(T) is slightly greater than A c [= A(T C )], we take only the first Fourier harmonics for the 
solution because the next harmonics are smaller near T c : 

u x (f) = u + 2uicos (-fi-T^j , (25) 
u y (f) = 2vi sin (jj^^ ■ (26) 



Exploiting the fact that w n 's are zero except for w and u± in Eq. (0), we get 

"° = 1 + ir + a(a + 2,^,0 ■ < 27) 

«i = - 2 W o- (28) 
Setting m = 2 in Eq. fl27|) and solving for A = A c , we have 

A c = Mai + ^e) (2Q) 

U 2 r]l(l + cie) 2 + 4af - 7T77i(l - de) 



Using the relations in Eqs. ( [T^D and , we plot A c against c<JoT r for different values of 
e in Fig. |T]. The maximal values of the crossover periods are more pronounced in the 
limit of smaller e, and A c 's converge to 2wai(= \/2ti) in the Hall regime and to 2n 2 r]i(= 
\f2nujQT r / y/e) in the dissipative regime. The reduced action integration near the crossover 
temperature can also be simply obtained by summing only n = and n = 1 contributions: 



J hd = 77A 
o 



9 / 27r 2 77i 47r 2 a? \ 9 
u 2 Q + 2 1 + — ii + — — i it 2 

u V A A ( A + 27r c i er ?i) / 



(30) 



which is reduced to 

/hd = ^A M 2 (3-m ), (31) 
o 

by using Eqs. (p7|) and (|28|) . 

The action integration Jhd obtained by the numerical variational method around T c and 
the one by Eq. (|31|) are compared in Fig. |^. The two curves in the figure perfectly join at 
the crossover period, which implies that our numerical method gives the correct solution. 
The dynamical magnetization relaxation rate Q is given by Q = Ti/Se- In real experiments 
Q is extracted from the magnetization M(t) = Mq[1 — Qln(t/to)} 0. From Eq. (|T2"D, we 
have Q(T)/Q = 2a/2// H d, where Q = (7rn s L c i? 2 e 5/2 ) _1 by taking Ci = c 2 = 1. Then, at 
the crossover temperature T c , since Jhd = 2A c /3, we have 

Q(T C ) _ SV2 

~qT~^v7' (32) 

In Fig. |3], Q(T C )/Q and Q(0)/Qo are plotted with respect to u r r . It is interesting that 



the shape of Q(T C )/Q for each e is close to that of Q(0)/Qo at zero temperature |l3| . This 



fact can be understood by considering the following features of the relaxation rate. Since 
the tunneling rate T is expected to be almost temperature independent for T < T c and 
T ~ exp(— U/ksT) for T > T c , the crossover temperature is approximately given by the 
relationship U/(k B T c ) ~ S E (T = 0)/K. And Q(0) = K/S E (T = 0) ~ (k B T c )/U = Q(T C ). 
Hence, we obtain the analytic expression for the relaxation rate at zero temperature using 
Eqs. (HD and © 



Q(0) ^ 3 y/(l + e) 2 + 4(^0^6 - (1 - e) 



by taking c\ — 1 and using ai/771 = lO^r/K^fe. Eq. ( p3|) agrees with the result in Ref. |^5j up 
to a numerical factor. Since the analytic form for Q(0)/Qo is known to be 5/6 in the limit 
of u T r — > 00, |L4] the correct prefactor of Eq. (|33|) is 5/12 instead of 3/(27r). Although 
Q(0)/Qo g° es to infinity as c^ot,. — > in Eq. (|33|) , it actually does not diverge in that limit, 
because, by including the inertia term not considered in this work, the approximate form 
of the classical action in the limit becomes S^/ft ~ L c y/m v Vo£,e 5 ^ 2 , which is independent 
of u>QT r , where m v is the inertia mass of a vortex. [jl],^,[^] In general, the mass term is 
relatively small in the Hall and dissipative regime and can usually be neglected. 

In a moderately clean regime, for small values of e each curve in Fig. |3] has a minimum 
around ujqt t = 1, which is interested. In fact, from Eq. ( |33|) we can see that the position 
of the minimum at zero temperature is u>oT r = (1 + e)/(l — e). As e becomes smaller, 
i.e., as j — > j c , the minimum becomes much more pronounced with its location moving 
toward ujQT r = 1 at the same time. The existence of such minima can be understood by 
considering the following qualitative features of the relaxation rates in the two regimes. Since 
u x (lo) is proportional to a\0j/ sinh(ai7rct;) in the Hall limit, Jl4| the classical trajectories with 
|co>| < 1/oti contribute to the Euclidean action mostly. From Eqs. ( [T2"D and (31) the correction 



to the Hall action S E H ^ by the small dissipation is given by (l+rji/a^S^ ~ (l+l/u> r r )S E , 
which leads to the relaxation rate given by Q(0) ~ ojqt t / (l + UoT r ). So the relaxation rate Q 
decreases with decrease in u T r from 00, which is also physically clear because the classical 
action increases by inclusion of the dissipation. In the opposite limit the correction to 
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the purely dissipative action Sg by the small Hall contribution is [1 + (ai/r]i) 2 ]S E ~ 
cuoTr[l + (&oT r ) 2 ], leading to the relaxation rate Q(0) ~ l/[a>oT r (l + (uor r ) 2 )]. In this case, 
Q(0) decreases with increase in u T r . Therefore, a minimum in Q(0) should exist in the 
intermediate regime, which suggests the existence of the strong pinning in the moderately 
clean samples. 



C. Quantum Relaxation in Dissipation Regime 

In the dissipative limit, we take «i — > in the action integration of Eq. (|24|) . The 
reduced action then becomes 

1 oo 

J hd = ^tA J2 0. + irrn\u n \)ul = I D . (34) 

n=— oo 

Noting that u n in the dissipative limit is given by u n = «oexp(— b\n\) where uq = 4ir 2 rji/A 
and b = tanh~ (27r 2 77i/A), pl| we get the reduced action given by 



' D 



1 / T * 21 



1 3 VT C 



(35) 



where A c = 2% 2 r]i and fcsT c = hy/cieVo / \2\f2n 2 R 2 a r]i) . In Fig. |L the relaxation rate using 
Eq. (|3li| ) is compared with the one obtained from the numerical solution: the two curves 
match quite well asymptotically in the region of small ujqt t values. 

D. Quantum Relaxation in Hall Regime 

In the Hall limit we take 771 — > in Eq. (0), which leads to the reduced action given by 

j 00 

J hd = -A £ (l + o?^)«2=/ H , (36) 



6 

n=— 00 



where u n satisfies 



00 

[I + a\u) 2 n )u n = " ^ U n+m U m . (37) 



n=— oo 
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While the instanton solution can be obtained analytically in the dissipative regime, the 
solution of Eq. fl3"T|) can be found numerically We use Eqs. ([16]) and (|T/]) rather than Eq. 
(|37|) and obtain the reduced differential equation 

d 2 u 

2aj-^ - 2u x + u 2 x = 0. (38) 

We then integrate for J H d with r\\ = in Eq. (|13|) using the solution for u x in Eq. (|38|) and 
obtain the reduced action In- As in the case of the dissipative regime, we have found that 
the relaxation rate agrees with the one obtained from the variational procedure in the limit 

Of U!qT t — > OO. 



IV. DISCUSSION AT FINITE TEMPERATURE 

We now consider the problem at the finite temperature in the whole regime, i.e., when 
the Hall and the dissipative dynamics are simultaneously present. The solutions for u x {r) 
and u y {r) in Eqs. ( |T8"D and ( p~9|) are obtained through w n 's and v n 's of Eqs. (^) and ( |2T|) 
which are numerically obtained by the variational method. In Fig. |5|, we show u x (f) and 
u y (f ) for various periods A which exhibit the typical trend of the classical trajectories as the 
period is successively shortened. The peak-to-valley amplitudes of u x {r ) and u y {f) decrease 
as the period gets shorter, eventually becomes flat at T c , i.e., u x (f) = 2 and u y (f) = 0. We 
subsequently calculate the reduced action ( |13"D via ( p4|) and the corresponding relaxation 
rate. The three-dimensional plot of Q(T)/Qq versus ujqt t and A for e = 0.01 is shown in 
Fig. p. We have also plotted Q(T)/Q against A for the different values of ojqt t in Fig. [7|. 
As can be seen in the figure, u>o T r = 1 is the boundary of the two different behaviors of the 
relaxation rate: for ujqTj. < 1, the relaxation rate increases with decreasing ioqt Ti whereas for 
ojQT r > 1 it increases with increasing u>oT r . The dependence of Q(T)/Q on u r r is shown in 
Fig. |8| for the different values of the period. 

What is interesting is the behavior of the relaxation rate in the intermediate region of 
u T r . We focus our attention on the four different temperature regimes, as indicated in Fig. 
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|8] (b). If the temperature is sufficiently low so that A > A( m ), the line yields no intersection 
points, and there exists quantum relaxation in the whole regimes of Uqt t . If the temperature 
is sufficiently high so that A < A^, quantum relaxation occurs only in the dissipative regime 
(ujQT r < (u)oT r )i)). In the temperature range A^ < A < A^ m \ on the other hand, quantum 
relaxation exists either in the Hall regime or in the dissipative regime, and purely thermal 
relaxation occurs in the crossover region between the two regimes. The values for , A W , 
and (u}QT r )[) can readily be computed from the position of the minimum and using Eq. ( |29| ) 
withci = 1: A( m > = ir(l+e)/y/2e, = ^tt, and {uoqT^d = y/e/(l-e). The corresponding 
relaxation rates are given by Q(A( m ))/Q = 6y/e/[ir(l + e)] and Q(A^)/Q = 3/vr. The 
minimum of Q(T)/Q in the intermediate regime is then noticed. As the temperature 
becomes lower, the quantum relaxation rate is more developed in the intermediate regime, 
and at an extremely low temperature it has a minimum at loqt t ~ 1. This feature is more 
pronounced for smaller e and larger A. Correspondingly, in such a regime the quantum 
depinning of a vortex is expected to be smaller at lower temperatures in the regime. 

Before concluding, we illustrate our results with specific numbers. In the experiment 
of Ref. 0, the relaxation rate is Q(0)/Q ~ 2.3 (2.0) in the YBCO (BiSCCO) system. 
In this case, ujqt t ~ 0.29 (0.37) for YBCO (BiSCCO) which corresponds to the Hall angle 
®h = arctan(a;o Tr) ~ 16° (20°), which depends on the oxygen content. The numbers imply 
that the samples are moderately clean. However, the regime which was considered in Ref. 
H was uj r r < 1, where the onset of the minimum just takes place. In order to observe the 
minima, the experiment should be extended to the region ujqt t 3> 1. 

V. CONCLUSIONS 

In conclusion, we have considered quantum tunneling of a vortex in the presence of the 
Hall and the dissipative dynamics. We have derived the analytic expression for the relaxation 
rate at zero temperature and obtained the numerical solutions by the variational method at 
finite temperatures. The relaxation rate is constant in the Hall limit and proportional to 
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l/(u T r ) in the dissipative limit, and, consequently, a minimum exist at u T r = 2{j c /j){\ + 



y 1 — 3 1 3c) — 1- Therefore, the strongest pinning is expected in the moderately clean sample 
at zero temperature. At finite temperatures, the quantum relaxation rate tends to vanish 
in the intermediate regime where both the Hall and the dissipative terms contribute to the 
dynamics of a vortex. At sufficiently low temperatures, quantum vortex tunneling occurs in 
the whole regime and the corresponding relaxation rate has a minimum at u T r ~ I. These 
features are expected to be observed in future experiments. 
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FIGURES 

FIG. 1. A c versus loqt t where c\ = 1. Note that the crossover temperatures become indepen- 
dent of e in the Hall regime {uqt t — * oo). 

FIG. 2. The relaxation rate Q(T)/Qo versus A near the crossover temperature for e = 0.1 and 
ujqt t =1, where A c ~ 7.6. The solid line represents the analytical curve from Eq. and the 

dotted line with diamonds indicates the result of the numerical calculation. 

FIG. 3. The relaxation rate evaluated at the crossover temperature: Q(T c )/Qq versus ujqt t 
where e = 0.1 (a), 0.01 (b), and 0.001 (c). Inset: the relaxation rate Q(0)/Qq at zero temperature 
with e = 0.1 (a), 0.01 (b), and 0.001 (c). 

FIG. 4. The relaxation rate in the dissipative limit when e = 0.1 and A = 10. The solid 
line represents the evaluation of Eq. (|35|), and the diamonds are the numerical results from the 
variational method. 

FIG. 5. Typical instanton solutions with different periods: u x (f) (top) and —iu y (f) (bottom) 
for e = 0.1 and ojqt t = 1.0, where the periods are oo (a), 10 (b), 8 (c), and 7.617 (d). 

FIG. 6. The relaxation rate Q(T)/Qq for e = 0.01 against UqT t and A. In order to show the 
curve for Q(T c )/Qq, we have omitted the purely thermal relaxation rate. See Fig. || for details. 

FIG. 7. The relaxation rate Q(T)/Qq versus A for the different values of uiQT r when e = 0.01. 
ujqt t increases from the bottom ( tvoT r = 1) and approaches oo. Inset: The case for ujqt t < 1. u>oT r 
decreases from the bottom. 
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FIG. 8. (a) Q(T)/Qq versus luqt t for different periods when e = 0.01. Periods are 30, 20, 15, 
and 10 (from the bottom). The curve with the period of 30 is already very close to that of an 
infinite period, which corresponds to zero temperature. The dotted curve is Q(T c )/Qq, of Fig. ||. 
In the region above the dotted curve, purely thermal relaxations exist along the horizontal lines, 
(b) A schematic diagram of A c versus ujqt t with the lines of constant temperatures. Four cases are 
considered : (I) A > K { c m \ (II) A = A^ m) , (III) A ( C H) < A < A^ m) , and (IV) A < A^ } . Note that 
Ai m) = 22.4, Q(T C )/Q = 0.189, and (uj T r ) D = 0.101. 
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